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Observations of Gamma-Ray Bursts with the Fermi Large Area Telescope have prompted theoretical advances 
and posed big challenges in the understanding of such extreme sources, despite the fact that GRB emission 
above 100 MeV is a fairly rare event. The first Fermi/LAT GRB catalog, published a year ago, presented 28 
detections out of 300 bursts detected by the Fermi Gamma-Ray Burst Monitor (GBM) within the LAT field of 
view. Building on the results from that work and on recent development in the understanding of the systematic 
errors on GBM localizations, we developed a new detection algorithm which increased the number of detections 
by 40%. Even more recently the development of the new event analysis for the LAT (’’Pass 8”) has increased 
the number of detections within the first 3 years of the mission to 45, up 50% with respect to the published 
catalog. The second LAT GRB catalog, in preparation, will cover more than 6 years of the mission and will 
break the barrier of 100 detected GRBs, a more than 20-fold improvement with respect to observations before 
the Fermi era in the same energy range. We will review the main features of the new algorithm, as well as 
preliminary results from this investigation. 


1. High-energy emission from 
Gamma-Ray Bursts 


Fermi/hAT observations are uncovering new and 
unexpected properties of the high-energy emission 
from Gamma-Ray Bursts (GRBs) shedding light on 
physics mechanisms, such as particle acceleration and 
emission processes, in u ltra-relativistic regime. The 
first LAT GRB catalog [Ackermann et alll2013 | con¬ 
tained 35 bursts, with 28 detected above 100 MeV 
with the standard likelihood analysis and 7 with the 
LAT Low-Energy technique (LLE). It established new 
high-energy features of GRBs, namely: 


• Additional power-law component during prompt 
emission: the prompt emission of most GRBs 
have been successfu lly described in th e past with 
the Band function [Band et al.lll993j| . In many 
bright GRBs observed by the LAT an additional 
power-law component is required to account for 
high-energy data. 


• Delayed onset: the emission above 100 MeV is 
systematically delayed with respect to the low- 
energy emission seen in the keV-MeV energy 
range. 


• Extended duration: the emission above 100 
MeV is also systematically longer than the 
prompt emission, and decays smoothly as a 
power law with typical decay index of —1, point¬ 
ing to a different physical origin with respect to 
the spiky prompt emission. 

Early afterglow models iKumar and Barniol DuranI 
|2009[ iGhisellini et al.l |2010| could explain the ob- 
served time decay, the delayed onset as the out¬ 
flow deceleration time scale, and the lack of variabil¬ 
ity. Hadronic models could explain these features as 


well: the onset delay could be the time to produce 


electromagnetic cascades [e.? 


Dermer and AtovanI 

l2006l iGuota and Zhan2ll2007l 

Asano and Inouel 2007| 


or the time required to accelerate, accumulate, 
and cool down relativistic protons via proton- 
synch rotron emission in a very strong magnetic 
field Razzaaue et al.l l2009j| . Eurthermore, energy- 
dependent delays would be expected in proton- 
synchrotron models where the cooling break shifts to 
lower energies at later times. On the other hand, 
we have analyzed the record-breaking GRB 130427A 
finding that, while several features of the “prompt” 
emission are in agreement with some internal-shocks 
models, other key details such as the estimated 
Lorentz factors for the colliding shells contradict 
exp ectations, possibly calling for different scenar¬ 
ios [Preece et al.l 12014 1. Also, the maximum energy 
of the photons in conjunction with a featureless high- 
energy light curve rules out both synchrotron and Syn¬ 
chrotron Self Compton from Fermi-accelerated parti¬ 
cles as emission mechanism, a radical departure from 
the s tandard external-shock model for the LAT emis¬ 
sion Ackermann et all 12014 1. New observations and 
insights are hence needed to foster further theoretical 
developments and overcome these difficulties. 

Our new specialized analysis described in Sect, [fal¬ 
lows the detection of faint high-energy GRB counter¬ 
parts, increasing the efficiency of detection by more 
than 50% and, when used in conjuction with the 
new “Pass 8” event selection, yielding more than 100 
bursts over the time span of the Fermi mission. The 
analysis and characterization of this new sample is in 
progress. When completed, it will provide the needed 
new insights. In particular, we will be able to settle 
some open questions which we could not firmly estab¬ 
lish in the first catalog due to the limited statistics: 


The existence of a separate population of hyper- 
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energetic events, characterized by a ratio be¬ 
tween high- and low-energy fluence much larger 
than the others. 


• The high-energy emission for all but 2 GRBs 
in the sample decays as t~^ at late times, 
which in the context of the fireball model fa¬ 
vors an adiabatic expansion. The other two 
GRBs instead deca .y as expected from a radia¬ 
tive regime (t~^ - ^) Blandford and McK^ll976L 


iKatz and PiranI Il997j . but they suffered from 
less-than-optimal observing conditions. Is there 
a class of truly radiative fireballs? 


• One of the GRBs in the sample presented a high- 
energy cutoff in the extra power-law component, 
while another one in the low-energy Band com¬ 
ponent. How common are spectral cutoffs? Do 
bursts with a cutoff have any other peculiar fea¬ 
ture? 


Analyzing the new enlarged sample will answer 
these questions, and will uncover new features. 


2. Pass 8 and GRBs 

The newly-developed “Pass 8” analysis is an event 
analysis rebuilt from the ground up with respect to 
previous versions (“Pass 6” and “Pass 7”). We refer 
the reader to other contributions in this Symposium 
for an introduction and for all the details about such 
analysis. Here we show how it will improve the num¬ 
ber of GRB detections. We started by performing a 
large number of simulations of GRBs. In particular, 
we fixed 10 flux values logarithmically spaced between 
10“® and 2 x 10“^ erg cm“^ s“^. Since the LAT ef¬ 
fective area is a function of the cosine of the off-axis 
angle 9 of the GRB, we also fixed a grid of 10 values 
uniformly distributed in cos 9. For each pair (flux, 
cos 9) we performed 200 simulations. These simula¬ 
tions were performed with the tool gtobssim, which 
accounts for the Poisson nature of the observation. 
For all simulations we used a fixed power-law spec¬ 
trum with index — 2, which is what is typ ically ob¬ 
served by the LAT [Ackermann et af1l2013j . We then 
overimposed these simulated GRBs on a background 
taken from a real observation, and performed a like¬ 
lihood analysis for each realization, recording the re¬ 
sults. We repeated this procedure using “Pass 7” data 
at first {P7REP_SOURCE_V15 class), and then a pre¬ 
liminary version of “Pass 8” data {P8.SOURCE_Vl 
class). We then compared the results. In order to 
determine the relative sensitivity between the two 
datasets, we computed the median significance for 
each point in the (flux, cos0) grid, computed as the 
square root of the TS value. The results are shown 
in Fig. [T] It is clear that Pass 8 can detect GRBs at 
much lower fluxes. We stress that the actual value 


for a particular detection threshold (for example the 
5cr one) depends on many factors (background level, 
duration of the GRBs, Zenith angle...) which are kept 
fixed in our simulations. Therefore, this is a compar¬ 
ative study, it is not a determination of the absolute 
sensitivity of Pass 8. 


3. New triggered search for Gamma-Ray 
Bursts 


Before the launch of Fermi it was estimated that 
the LAT would have observed 10-1 2 GRBs/year abov e 
100 MeV and 6-8 above 1 GeV [Band et fol l2009l| . 
During the first 3 years, however, observations were 
slightly below such exp ectations [Ackermann et all 
I 2 OI 3 L iGuetta et al.l 12011 1. Recently we developed a 
new search algorithm for GRBs, which is now up and 
running 24/7. As shown in Fig. [51 it provides 45% 
more detections than the algorithm used for the cata¬ 
log when using Pass 7 data, and 60% more when using 
the newly-developed Pass 8 data. This improvement 
was achieved exploiting the results of the first LAT 
GRB catalog, presented in Sect. [TJ The new algo¬ 
rithm consists of 10 searches running in parallel over 
time intervals logarithmically spaced from the trigger 
time to 10 ks after that. For each of these time inter¬ 
vals, these are the steps of the new algorithm: 


1. A trigger is received either in real-time through 
the GGN system^ or during an off-line analysis. 
In real-time most of the triggers are from the 
GBM, although also triggers from Swift, INTE¬ 
GRAL and other observatories are received and 
processed. 


2. If the trigger comes from the GBM, a “finding 
map” is produced to account for the position un¬ 
certainty. It has been recently reported that the 
GBM loca lizes GRBs with a systern atic error of 
up to 15° Connaughton et aL[l2014[| . Since this 
systematic error dominates over the statistical 
one, we always use finding maps of 30 x 30 deg. 
The finding map is essentially a TS map like the 
one produced by the tool gttsmap. Indeed, a grid 
in equatorial coordinates and with a spacing of 
0.7 deg is created covering the finding map. For 
each point of the grid, a likelihood analysis is 
performed including a new source at that po¬ 
sition. However, points in the grid having less 
than 3 photons within 10 deg are not considered. 
This avoids running a likelihood analysis, which 
is computer-intensive, in regions of the finding 


^ Gamma-Ray Coordinate Network, 

http://gen.gsfc.nasa.gov/ 
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Figure 1: Median significance as a function of flux and off-axis angle for Pass 8 (left) and Pass 7 (right). The dashed 
line is the 5cr threshold. Pass 8 goes to much lower fluxes than Pass 7. 


map where there is no possibility of detecting a 
new point source. At the end of each likelihood a 
TS value is computed, and associated with the 
point in the grid. When all points have been 
processed, the maximum of the TS in the map 
is considered the best guess for the position of 
the new transient, and marked for further anal¬ 
ysis. Note that if the trigger comes from Swift 
or INTEGRAL, this step is not executed, since 
their localization errors are much smaller than 
the typical size of the LAT PSF. 

3. The position of the candidate transient is opti¬ 
mized with the tool gtfindsrc 

4. A new likelihood analysis is performed on the 
best position found. If the TS from this final 
analysis is above a certain threshold, we consider 
it a new detection 


5. If running in real time, the results from the anal¬ 
ysis are used by Burst Advocates to disseminate 
alerts to the community through GCNs. 


In all likelihood analyses involved in the sequence 
the likelihood model consist of the Galactic template 
and the Isotropic template provided by the Fermi 
collaboration, plus all point sources from the 2FGL 
source catalog. While the normalization of the two 
templates are left free to vary, all parameters for the 
2FGL sources are kept fixed. If the algorithm finds 
a new candidate source with a position compatible 
with one of the 2FGL sources, further analysis are per¬ 
formed manually to distinguish between a real tran¬ 
sients, and other phenomena such as AGN flares. The 
algorithm does indeed trigger on flaring sources, which 
are al so found by othe r real-time a lgorithms such as 
ASP [Band et al.ll2n09l | and FAVA Ackermann et al.] 
l2013f . 


The procedure involves a certain number of trials, 
which might appear large. However, all the time win¬ 
dows are overlapping and involve the same region on 
the sky, therefore these trials are not independent 
at all, with many photons present in many of the 
time windows. Thus, the effective number of trials 
is rather small. While the full characterization is still 
in progress, preliminary simulations indicates that a 
5(7 detection corresponds to TS ~ 28. 

A first run of the new algorithm on ~ 5 years 
of preliminary Pass 8 data returned 86 detections. 
Adding the ^ 20 LLE-only detections, we have al¬ 
ready reached the milestone of 100 LAT-detected 
GRBs. Given the current rate of detections we also 
expect to exceed 100 likelihood detections within the 
year, an impressive milestone significantly exceeding 
pre-launch expectations. As expected, the new detec¬ 
tions populate the lower part of the fluence distribu¬ 
tion, demonstrating that the new algorithm and “Pass 
8 ” improve the sensitivity of the search. This is shown 
in Fig. [3l 

We also note that this preliminary study is using the 
“Pass 8” Source class for all time scales. It has been 
shown in the past that for short time scales (< 200s) 
using classes with larger acceptance at the expense of 
a larger background is beneficial. We will be using 
such classes in the final study, further increasing the 
number of detections. 
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Figure 2: Cumulative detections for the time span covered by the LAT GRB catalog [Ackermann et al]l2013l ]. The new 
analysis yields 45% more detections than the one used in the catalog when run on the same data (blue and green), and 
60% more with Pass 8 data. 



Fluence (10 keV -1 MeV) [erg cm“^ ] 

Figure 3: Left panel: low-energy fluence distribution for the whole GBM sample (upper panel) and for the sample 
detected with the new a lgorithm in LAT data (lower panel). Right panel: similar plot from the first LAT GRB catalog 
[Ackermann et aDl2013l| . The blue dashed line marks the region where most of the new detections lie, showing the 
increased sensitivity of the new algorithm with respect to the old one. 
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